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MASSACHUSETTS 


ABSTRACT 


The  large  proper  motions  and  diurnal  parallaxes  of  earth-bound 
artificial  satellites  coupled  with  the  use  of  devices  more  sophisticated 
than  the  photographic  plate  to  record  their  motion  force  one  to 
consider  methods  of  analysis  beyond  the  traditional  ones.  In  particular, 
it  is  theoretically  possible  to  deduce  from  the  streak  made  by  the 
passage  of  a  celestial  object  a  complete  specification  of  its  location 
and  velocity.  Thus,  for  artificial  satellites,  one  has  the  possibility 
of  orbital  element  set  construction  from  a  complete  set  of  initial 
conditions  without  recourse  to  any  approximations.  In  this  report,  I 
show  how  and  why  this  is  possible.  The  practical  situation  is  less 
hopeful  however,  in  that  the  distance  and  radial  velocity  can  be 
determined  only  if  the  streak  has  measureable  curvature.  Extremely 
accurate  (+  0'.'05/sec)  angular  velocities  should  be  obtainable.  The 
positional  accuracy  is  a  function  of  the  driving  of  the  telescope. 
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I.  MOTIVATION 


If  one  uses  a  solid  state  device  camera  (CID  or  CCD)  or  an  electron 


beam  tube  with  digitized  output  to  observe  artificial  satellites,  then  a 


whole  new  field  of  optical  observing  and  data  reduction  is  thrown  open 


With  such  a  device  connected  to  a  telescope  moving  at  the  sidereal  rates 


the  time  history  of  the  motion  of  a  satellite  can  be  obtained.  That  is,  we 


can  know  the  coordinates  of  the  light  deposited  on  the  camera  target  and 


when  that  light  was  deposited.  This  is  possible  with  a  photographic  plate 


but  difficult  to  do  rapidly  and  accurately.  Within  this  knowledge  is 


a  set  of  complete  information  concerning  the  satellite's  geocentric  location 


and  velocity.  Hence,  an  analysis  of  the  image's  location  as  a  function  of 


time  can  yield  a  full  set  of  initial  conditions  for  orbital  analysis 


Moreover,  there  should  be  an  extremely  high  internal  accuracy  in  these  points 


as  they  all  suffer  from  nearly  the  same  systematic  errors.  The  random  and 


systematic  errors  associated  with  the  telescope  pointing  enter  weakly  or 


not  at  all.  Hence,  we  have  the  promise  of  more  information  than  just  the 


object's  position  and  of  extremely  high  accuracy  in  this  additional  data 


The  analysis  of  the  imaged  motion  is  straightforward.  It  is  given 


in  Sections  II,  III,  and  IV.  In  Section  V,  a  practical  discussion  of  the 


situation  reveals  that  the  distance  and  radial  velocity  can  only  be 


determined  if  the  motion  on  the  camera  target  deviates  from  a  straight 


line.  The  angular  velocity  will  be  very  precisely  ascertained  and  this 


Information  should  be  used  in  initial  orbit  construction  as  well  as  in 


the  differential  correction  of  orbits 


II.  THE  GEOCENTRIC  GEOMETRY 

We  use  Che  standard,  almost  Inertial,  spherical,  geocentric  coordinate 
system  based  on  the  earth's  equator  and  the  Vernal  Equinox.  A  celestial 
object's  location  is  given  by  r  «  r_£,  Z  ■  (cos6cosa,cos<5sina,sin6) . 

The  geocentric  distance  is  r,  the  geocentric  declination  is  5,  and  the 
geocentric  right  ascension  is  a.  Over  very  short  times,  very  much  less 
than  an  orbital  period  (<0.3%  in  practice),  we  regard  the  motion  as 
uniform  and  rectilinear  in  this  coordinate  system.  Since  this  is 
different  from  uniform  "rectilinear"  (e.g.,  along  a  great  circle)  motion 
on  the  celestial  sphere,  there  is  a  coupling  between  the  quotient  of 
radial  velocity  (r)  and  the  distance  (r)  and  the  object's  position 
(e.g.,  Z).  To  see  how  this  comes  about  consider  the  orthonormal  basis 
given  by  l,  a,  and  d;  a  *  (-sina,cosa,0) ,  d  -  (-sin<5cosa,-sin6sina,cos6) . 
We  can  write 

i  ■  Z»r,  w  -  Acosfia  +  6 d.  (1) 

In  terms  of  the  constant  unit  vector  n  ■  (0,0,1),  which  points  towards 
the  North  Celestial  Pole, 

f  -  f.r,  (2a) 

&cos6  ”  a* 0)  ■  «•(£  x  r)/(rcos6),  (2b) 

6  ”  ”  (w  .fsecfi  -  ftan6)/r.  (2c) 
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By  differentiating  Eqs.  (2)  and  enforcing  the  constraint  V  =  £,  we  find. 


ru) 


(3a) 


a  =  2a(6tan6  -  r/r). 


(3b) 


6  =  -  &^sin6cos6  -  26(r/r), 


(3c) 


Thus,  even  though  the  space  motion  is  along  a  geodesic,  the  position  of  the 
object  as  a  function  of  time  changes  because  of  its  purely  radial  motion. 
This  term  is  known  as  the  foreshortening  term  in  astrometry.  It  can  be 
measured  for  some  natural  celestial  objects.  To  convince  the  reader  that 
r/r  is  not  numerically  much  smaller  than  &  or  6  for  artificial  satellites 
consider  a  satellite  with  argument  of  perigee  =  270°,  inclination  =  60°, 
eccentricity  *  1//2.  These  numbers  are  typical  of  the  Molniya  class 
satellite  of  which,  with  their  associated  rocket  bodies,  there  are  hundreds. 
For  such  a  satellite  observed  at  the  equator  4  =  n^2,  6  =  n/6 ,  r/r  =  2n 
where  n  is  the  mean  motion  (for  which  I'll  use  2  rev/day  when  it's  needed). 
Clearly,  the  three  quantities  are  comparable  and,  since  6=0,  the  fore¬ 
shortening  terms  drive  the  acceleration  of  the  position. 

To  derive  the  rigorous  results  we  use  the  equations  of  motion 


r (t)  =  r(t  )  +  r(t  ) (t  -  t  ), 
—  —  o  —  o  o 


(4a) 


K t)  -  t( tQ). 


(4b) 
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By  manipulating  the  component  forms  of  Eqs.  (A),  we  can  obtain 


(d  -  r(t)/r(t  ),  y  =  &(t  ),  y,  =  6(t  ),  v  =  r(t  )/r(t  ),  T  =  t  -  t  ,  6 
o  a  ooo  oo  oo 

=  6(t  ),  a  =  a(t),  etc.)* 
o 

d2  =  1  +  2vT  +  (v2  +  uj2)T2,  (5a) 


tan (a  -  aQ)  =  yaT/[l  +  (v  -  y6tan6Q)T], 


y*T  +  (1  -  d  +  vT)  tan<5 

tan{  (<5  -  6  )/2]  =  — - 2—r 

d  +  {[1  +  (v  -  y^tanfi  )T] 

0  o 


r/r  =  d/d  =  [v  +  (v2  +  U)2)T]/d2, 


2_2-i  1/2  ’ 

+  yaT  }  ' 


(5b) 


(5c) 


(6a) 


d  =  Ua/{U  +  (v  -  y6tan6Q)T]2  +  yV}, 

.  ^(1  +  vT)  "  w2Ttan6Q 

d2Ul  +  (v  -  y6tan6o)T]2  +  yV}1/2‘ 


(6b) 


(6c) 
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These  equations  were  first  correctly  derived  by  Eichhom  and  Rust  . 

Equations  (5  and  6)  give  a  complete  description  of  the  motion  in 
the  form  we’ll  need  it.  The  position  at  any  time  depends  on  the  five 
initial  conditions  aQ,  6q,  v,  y^,  and  y^.  To  recover  the  sixth  quantity 
we  need  to  refer  the  motion  to  the  observer's  location  instead  of  the 
center  of  the  earth.  Before  we  do  this,  we  could  put  the  physics  back 
into  the  problem.  My  main  purpose  in  considering  uniform  rectilinear 
motion  was  to  bring  to  the  reader's  attention  the  geometrical  effect 


A 


* 


J 


that  the  foreshortening  terms  do  influence  the  position.  Equations  (5  and  6) 
could  now  be  replaced  by  the  formulas  of  two-body  Keplerian  motion.  The 
development  of  the  next  two  sections  is  independent  of  the  physics,  or  lack 

thereof,  in  Eqs.  (5  and  6).  It  should  also  be  clear  that,  since  r/r  enters 

•  •  *• 

first  into  a  and  6,  that  it  can  only  be  determined  if  the  imaged  path  has 
measureable  curvature.  We'll  calculate  the  radius  of  curvature  in  §V  and 
find  it  to  be  essentially  infinite.  Hence,  as  a  practical  matter,  the 
physics  doesn't  matter. 


III.  DIURNAL  PARALLAX 

I  shall,  except  for  time  variables,  uniformly  use  the  corresponding 

upper  case  letter  to  denote  a  topocentric  quantity.  Thus,  A  =  A(t)  = 

the  topocentric  right  ascension  at  time  t,  =  Aq  =  A(tQ)  =  the  proper 

motion  in  the  topocentric  declination  at  time  t  ,  D  =  R/R  etc.  Let 

o  o 

£=  p(cos<}>  'cost,cos<J>'  sinx,  sin<J>' )  be  the  observer's  geocentric  location. 
Here  p  is  the  observer's  geocentric  distance,  <p'  his  geocentric  latitude, 
and  T  is  the  local  mean  sidereal  time.  The  translation  to  the  observer 
is  accomplished  by 


£  =  R  +  £. 


(7) 


If  we  write  the  right-hand  side  of  this  in  component  for,  and  use  Eqs. 
(4)  similarly  expressed,  one  can  derive  the  results  analogous  to  those 
in  Eqs.  (5  and  6).  To  persuade  the  reader  that  we  don't  want  to  follow 
this  course  one  finds,  for  the  location, 


D2  +  2DScosZ  =  1  +  2VT  +  (V2  +  fi2)T2  +  2ScosZ 

o 


+  2ST(VcosZQ  +  dcosZQ/dt)  +  (StTcos<f> ' ) 2 


-  2StT2cosif)'cosA  sinH  (V  -  MxtanA  -  M  cotH  ) ,  (8a) 

o  o  o  o  a  0 


M  T  +  S 

tan  (A  -  A  )  =  - ^ - - - ,  (8b) 

1  +  (V  -  M.tanA  )T  +  S„ 

0  0  2 
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tan[(A  -  A  )/2]  =  [M.T  +  (1  -  D  +  VT)tanA  ]/(D  +  C) .  (8c) 

O  0  o 


The  missing  initial  condition  is  S  =  p/R^.  The  other  quantities  are 
defined  by  (Z  =  topocentric  zenith  distance,  H  =  topocentric  hour  angle) 


cosZ  =  sin^'sinA  +  costjj'cosAcosH, 


(9a) 


H  =  t  -  A, 


(9b) 


C2  =  [1  +  (V  -  M-tanA  )T]2  +  M2T2  +  2{S,M  T 
6  o  a  la 


+  S2[l  +  (V  -  M^tanAQ)T] }  +  (Scos<j>' secAQ) 2{  (tT) 2 


+  2 [ 1  -  cos(t  -  t  )  -  fTsin(x  -  T  )]}, 
o  o 


(9c) 


S,  =  Scos<b'secA  IsinH  +  tTcosH  -  sin(T  -  A  )], 
1  o  o  o  o 


(9d) 


S_  =  Scoscb'secA  [cosH  -  tTsinH  -  cos(x  -  A  )].  (9e) 

2  o  o  o  o 


Necessarily  as  S  +  0  D  +  d  etc.  and  Eqs.  (8)  reduce  to  Eqs.  (5).  The  analog 
of  Eqs.  (6)  can  be  obtained  from  Eqs.  (8)  by  differentiation  with  respect 
to  the  time.  Equations  (8  and  9)  are  rigorous. 

It  seems  simpler  to  me  to  interpret  Eqs.  (7)  geometrically.  Then  one 
f inds^ 
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■ 


tan  (a  -  A)  =  tan(H  -  h)  =  - ^cosf  sec6/d)  sinh  .  ,  (10a) 

II  -  (scos<J>'sec6/d)cosh] 


tan(6  -  A) 


(ssintft1 cscy/d) sin(y  ~  6) 

II  -  (ssin<t>'cscY/d)cos(Y  -  6)]’ 


(10b) 


D  =  dsin(6  -  y)csc(A  -  y) ,  (10c) 

where 

tany  =  tan<J> '  cos  I  (a  -  A)/2]secIh  +  (a  -  A)/2].  (lOd) 


Equations  (5,  10)  completely  specify  the  topocentric  location  of  the  object 

in  terms  of  the  initial  conditions  s  =  p/r  ,  a  ,  6  ,  v  =  r  /r  ,  y  =a, 

o’  o’  o’  o  o’  Ha  o’ 

and  =  6q.  When  the  telescope/camera  combination  images  the  motion,  with 
the  telescope  in  sidereal  drive,  on  a  photographic  plate  or  camera  target, 
the  resulting  streak  contains,  albeit  implicitly,  a  complete  description 
of  the  motion.  The  last  step  then  is  to  project  a  portion  of  the 
celestial  sphere  onto  a  plane. 


IV.  STANDARD  COORDINATES 

If  the  camera  is  equivalent  to  a  pinhole  camera,  then  the  projection 
of  the  imaged  portion  of  the  celestial  sphere  onto  a  plane  will  be  a 
gnomonic  projection.  The  resulting  coordinate  system  on  the  plane  is  a 
rectangular  one  known  as  the  standard  coordinate  system.  To  relate 
standard  coordinates,  symbolized  by  £,  n,  to  right  ascension  and  declina¬ 
tion  we  choose  the  unit  of  length  for  the  standard  coordinates  to  be  the 
telescope's  focal  length  and  we  need  to  know  the  coordinates  where  a 
plane,  parallel  to  the  plane  of  the  plate,  touches  the  celestial  sphere. 
This  is  known  as  the  tangential  point.  One  end  of  the  optical  axis  of 
the  telescope  pierces  the  celestial  sphere  at  this  point,  the  other  end 
intersects  the  center  of  the  plate.  The  positive  n  axis  points  toward 
the  North  Celestial  Pole.  The  positive  £  axis  points  toward  the  east 
point^. 

If  (a*,  6*)  are  the  coordinates  of  the  tangential  point  then  the 
relationship  between  £,  n  and  a  corresponding  a,  6  is  given  by 

£  =  sec6*sin(a  -  a*) /[tan6*tan6  +  cos(a  -  a*)],  (11a) 

n  =  Itan6  -  tan6*cos(a  -  a*) ] /[tan6*tan6  +  cos (a  -  a*)].  (lib) 

The  inverse  of  Eqs.  (11)  is 

cot6sin(a  -  a*)  "  ^/(sinS*  +  ricosfi*),  (12a) 
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cot6cos(a  -  a*)  =  (cos6*  -  nsin6*)/(sinS*  +  ncosd*) .  (12b) 

The  formula, 

sin6  ■  (sin6*  +  ncos<$*)/(l  +  Z?  +  (12c) 

is  useful  near  the  equator. 

The  problem  is  now  completely  solved.  Since  we  obviously  observe 
topocentrically,  we  replace  a,  6  in  Eqs.  (11)  by  A,  A  computed  from 
Eqs.  (10).  The  geocentric  coordinates  appearing  in  Eqs.  (10)  are  given 
as  a  function  of  time  by  Eqs.  (5).  Since  we  measure  the  properties  of 
the  streak  imaged  on  the  plate,  we  have  an  implicit  problem  for  the 
determination  of  the  location  and  velocity  of  the  moving  object.  Once 
these  have  been  solved  for  we  can  compute  the  orbital  elements.  It's 
also  clear  from  Eqs.  (11,  12)  that  it  is  simplest  if  the  position  of 
the  telescope  is  not  changed  during  the  exposure  of  the  plate. 
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V.  APPROXIMATE  FORMULA 

By  combining  all  of  the  formulas  of  the  preceding  three  sections,  we 
find  that,  through  all  terms  of  the  second  order 

K  =  {Ea  +  MaT  +  [2Ma(M6tanAo  -  V)  +  S^cosij)’  secAosinHo]T2/2}cosAo 

-  M6T(Ea  +  M6T)sinAQ,  (13a) 

n  =  E.  +  McT  -  [M2sinA  cosA  +  2VM.  +  Si2cosd>' sinA  cosH  ]T2/2 
oo  aoo  0  oo 

+  [(Ea  +  MaT)2/2]sinAocosAo,  (13b) 

where 


E  =  A  -  a*,  Er  =  A  -  6*. 
a  o  o  o 


(13c) 


The  reader  is  reminded  again  that  these  expressions  are  not  series  in  S 

but  rather  series  in  Ea>  Eg,  M^T,  MgT,  and  VT.  It's  clear  that  all  of 

these  quantities  are  limited  by  the  f ield-of-view  of  the  telescope/camera 

combination.  It  should  also  be  clear  that  neither  a  -  a*  nor  6  -  6* 

o  o 

is  limited  by  the  field  of  view  and  that  they  aren't  necessarily  small. 
(Compare  with  references  5,  6,  and  7.) 

The  first  thing  we  shall  compute  is  the  curvature  of  the  streak.  If 
the  radius  of  curvature  is  effectively  infinite,  then  the  quadratic  terms 
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in  the  time  can  be  dropped  from  Eqs.  (13)  and,  as  a  practical  matter,  the 
distance  and  radial  velocity  are  unobtainable.  The  leading  term  in  the 
curvature  is 

K  ■  (St^cos^'/J)^)  iM.sinH  +  M  sinA  cosA  cosH  |.  (14a) 

1  6  o  a  o  o  o 

As  S  +  0  the  curvature  becomes 

K  -*•  2 1  ( VEa/«)  sinAQ  | .  (14b) 

For  the  satellite  considered  earlier,  if  its  equator  crossing  occurs 

on  the  meridian,  K  'v  1/400.  For  any  reasonable  focal  length  the  streak 

will,  therefore,  be  indistinguishable  from  a  straight  line. 

The  high  internal  accuracy  of  the  measured  values  of  n  imply 

that  E  cosA  ,  Ex,  M  cosA  ,  and  M.  will  be  very  well  determined.  The 
oi  o  o  ol  o  o 

external  systematic  errors  in  the  angular  velocity  will  also  be  negligible. 
However,  the  external  systematic  and  random  errors  in  the  position  will 
reflect  those  of  the  telescope  pointing.  In  fact,  in  the  simplest  situa¬ 
tion,  the  standard  deviation  of  E  cosA  and  Ex  will  be  /N  down  from  those 

a  o  6 

of  £  and  n  (W  =  number  of  observations)  while  the  standard  deviations  of 
MacosAo  and  will  %  (total  exposure  duration)  down  from  these.  This  may 
well  be  a  few  hundredth's  of  a  second  of  arc/second.  The  length  of  the 
streak,  when  K  *  0,  is  just  ft  times  the  duration  of  the  exposure.  This 
may  be  useful  as  a  constraint  but  because  of  its  relatively  large 
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standard  deviation  is  not  useful  for  the  overall  analysis.  Finally, 
had  only  the  developed  streak  been  analyzed,  we  can  see  from  Eqs.  (13), 
that  the  model  would  have  been 

£  -  (EaMfi  -  EfiMacosAo)/M6  +  O^cosA^M^n.  (15) 

Hence,  unless  the  streak  passes  through  the  origin  and  the  length  of  the 
streak  is  employed,  all  of  the  Information  would  have  been  lost. 
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ABSTRACT  (Cornu mum  on  revert,  tid.  if  nectary  and  identify  by  block  number) 

The  large  proper  motions  and  diurnal  parallaxes  of  earth-bound  artificial  satellites  coupled  with 
the  use  of  devices  more  sophisticated  than  the  photographic  plate  to  record  their  motion  force  one  to 
consider  methods  of  analysis  beyond  the  traditional  ones.  In  particular,  it  Is  theoretically  possible 
to  deduce  from  the  streak  made  by  the  passage  of  a  celestial  object  a  complete  specification  of  Its 
location  and  velocity.  Thus,  for  artificial  satellites,  one  has  the  possibility  of  orbital  element  set 
construction  from  a  complete  set  of  initial  conditions  without  recourse  to  any  approximations.  In 
this  report,  I  show  how  and  why  this  Is  possible.  The  practical  situation  Is  less  hopeful  however,  in 
that  the  distance  and  radial  velocity  can  be  determined  only  If  the  streak  has  measure  able  curvature. 
Extremely  accurate  (gr.'oS/seft)  angular  velocities  should  be  obtainable.  The  positional  accuracy  is  a 
function  of  the  driving  bf  the  telescope. 
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